HYDRODYNAMIC LIMIT OF SYMMETRIC EXCLUSION PROCESSES 

IN INHOMOGENEOUS MEDIA 



A. FAGGIONATO 



Abstract. In [j] M. Jara has presented a method, reducing the proof of the hydrody- 
namic hmit of symmetric exclusion processes to an homogenization problem, as unified 
approach to recent works on the field as [N] , |F1| , |F2) and [FJL] . Although not stated in 
[j] , the reduction of the hydrodynamic limit to an homogenization problem was already 
obtained (in a different way) in [N], |F1) . This alternative and very simple relation be- 
tween the two problems goes back to an idea of K. Nagy [N], is stated in [Fl] [Section B] 
for exclusion processes on if" and, as stressed in [F2] . is completely general. The above 
relation has been applied to [N], [Fl] . [F2) and [FJL] and could be applied to other 
symmetric exclusion processes, mentioned in [j]. In this short note we briefly recall this 
unified approach in a complete general setting. Finally, we recall how the homogenization 
problem has been solved in the above previous works. 
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1. Introduction 

Exclusion processes in inhomogeneous media can be used to analyze transport properties 
of particle systems in inhomogeneous and/or disordered media with hard-core interaction. 
Due to the possible lack of spacial invariance, the resulting model is an example of non- 
gradient exclusion process, in the sense that the transition rates cannot be written as 
gradient of some local observable. Despite this fact, if the jump rates do not depend 
on the orientation of the jump (as in [Nj, [FT], |F2j . [FJLj but not as in Ql , |Q2| and 



|FM| ) ■ the hydrodynamic limit of the exclusion process can be proven without using the 
very sophisticated techniques developed for non-gradient systems (cf. [KL] and references 
therein), which in addition would require non trivial spectral gap estimates that fail in 
the case of jump rates non bounded from below by a positive constant (cf. Section 1.5 in 
[M]). The strong simplification comes from the fact that the infinitesimal variation of the 
occupancy number at a give site is a linear combination of occupancy numbers. 

Alternative routes to prove the hydrodynamic limit are the method of corrected empir- 
ical measures developed in [JLj and |GJj . which is also at the base of fJ], and a duality 
relation developed in [N], |Flj . The advantage of the former is that it works also for 
zero range processes, while the latter is very elementary and direct. In both cases, the 
hydrodynamic limit reduces to an homogenization result for a random walk (possibly in a 
random environment). In [J] the connection between hydrodynamics and homogenization 
by means of the corrected empirical measure has been described in a general setting, and 
presented as "unified approach" to recent works on the field as [N], |Flj . |F2j and |F JLj . 
The alternative connection due to the duality relation presented in [N] and |Flj is not 
mentioned at all in [T. On the other hand, this is indeed the connection between hy- 
drodynamics and homogenization used in [Nj, [Flj . |F2j and [FJLj . and could be applied 
as well to the other symmetric exclusion processes. In what follows we briefly recall this 
equivalence working in a more general setting than [Nj, [Flj . 
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Finally, we remark that, having established the equivalence between the two problems, 
it remains to solve the homogenization problem. In [N] the author proves indeed an 
invariance principle for the random walk, under rather restrictive conditions. In |F1| the 
author applies the method of Stone [S] for random walks that can be written as space- 
time change of Brownian motion, solving the hydrodynamic limit for symmetric exclusion 
processes on Z under very natural conditions. Again, Stone's theory has been applied in 
|F JL| . In [F2| the homogenization problem has been solved by means of the notion of 
two-scales convergence ( |Nu] . [S] . [ZP] ) and is indeed the core of |F2j . In particular, in 
[j] the author does not recover the previous results of [N], [FT], |F2) and [FJLj as stated 
in the Introduction of [Jj, since the homogenization problem is not solved in [j]. 

2. Models and main result 

For each n ^ 1 let Gn = (Vn, En) be a graph with unoriented bonds. Vn and En denote 
respectively the family of vertexes and the family of edges. We suppose that Vn is a locally 
finite subset of X, where X is a metric space endowed of a Radon measure m (one could 
as well simply require that Vn can be injected in X with locally finite image, this would 
slightly change the notation below). Moreover, we assume that 

rUn := — 6x -i- m (2.1) 

with respect to the vague topology, for a fixed sequence a„ such that 

lim ttn = oo . (2-2) 

n—>-oo 

The state space of the exclusion process on G„ is given by the set {0, 1}^", endowed of 
the product topology. 

The exclusion process on the graph Gn is defined in terms of a rescaling constant 6„ and 
a conductance field {c„(6) : b G En}- It is a Markov process on {0, 1}^" whose Markov 
generator acts on local functions (i.e. depending only on ijx with x varying in a fixed finite 
subset of Vn) as 

Cnf{v) = bn Yl c„(x,y)(/(7?^'^)-/(7?)). (2.3) 

{x,y}&En 

Above we have written c{x, y) for the conductivity c({x, y}) associated to the bond {x, y} 
and we have used the standard notation 

r]y \i z = X , 
rj^ if z = y, 
r]z otherwise . 

We assume that the above dynamics is well defined (cf. [L] for sufficient conditions). It 
can be constructed as follows. On a common probability space (il, J-", P) define a family 
of independent Poisson processes NJ^i-), parameterized by 6 G En, such that E{N^{t)) = 
Cn{b)t. Roughly, if t is a jump time of A^^ then at time t the exclusion process performs 
the jump r]{t—) i— )• rj{t)^. To formalize the above definition we assume that there exists 
e > such that for almost all uj the connected components of the graph ^"(w), obtained 
from Gn by keeping only the bonds b such that N^{e)[uj] ^ 1, have finite cardinality (cf. 

In what follows, given t ^ we denote by ^"(w) the random graph with vertexes Vn 
and edges b e En such that N^{t)[uj] > Nnike)[u;], where ke ^ t < {k + l)e, A; G N. 
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As consequence of the above assumption, there exists a measurable subset A & with 
P(^) = 1, such that Qt{ijj) has only connected components with finite cardinality for 
all t ^ and all u ^ A. Let u ^ A. Then, given an initial configuration 7y(0), the 
configuration 'q{t) = 'q{t)[uj] at time t is defined as follows: 
Let C be any connected component of Qt'{ijo) and let 

{si < S2 < • • • < Sr} = {s : N^{s) = Kis-) + 1, b bond in C, < s ^ t}. 

Start with r/(0). At time si switch the values between rj^ and r/y if 6 = {x,y}, NJ^{si) = 
N^(si—)+l and 6 is a bond in C. Repeat the same operation orderly for times S2, S3, ... , Sr- 
Then the resulting configuration coincides with on C. 

Below, we denote by P^^ the law of the exclusion process on G„ when the initial 
distribution is 

As special case one can consider the dynamics of a single particle. The resulting process 
is a continuous-time random walk on G„, having Vn as state space and infinitesimal 
generator 

hnh{x) = bn Yl Cn{x,y){h{y)-h{x)) . (2.4) 

Note that, given vertexes a,b € Vn, we write a ~ 6 if {a, 6} G En- We write Xn{t\x) for 
the above random walk when starting at site x. 
Given a function (p on Vn, we write 

Pr^{x):=E[ip{Xn{t\x))], xeVn. 

Theorem 2.1. Suppose that there exists a stochastic process W on X, with symmetric 
probability kernel, satisfying the following homogenization property: 

lim — V IPr^pix) - PMx)\ = , V99 G C,{X) , (2.5) 

where 

PMx) :=E[^{W{t\x) )] 

and W{t\x) denotes the process W starting at x € X. Then the following hydrodynamic 
behavior of the exclusion process holds: 

Let po : X ^ [0, 00) be a bounded Borel function and let {pn)n ^ i be a sequence of 
probability measures on {0, 1}^" such that, for any 6 > and any G {Pt^^ : t Ojtp G 
Cc{X)}, it holds 



lim pn 

ntoo 



— ip(x)ri(x) — [ ip(u)po(u)m(du) > d] =0. (2.6) 



Then, for any 6 > 0, t > and and any ip G Cc{X) it holds 

1 



nfoo 



— f{x)v?{x) - I (p{u){Ptpo){u)m{du] 



>6]=0. (2.7) 



We point out that, in applications, the set {Pt'i' : t 0,ip & Cc{X)} is given by enough 
regular functions decaying to at infinity. 

One can restate the homogenization property (j2.5p in a form more related to partial 
differential equations as done in [F2j , by means of the concept of 2-scale convergence which 
is particular suited for the above setting (cf. [ZPj ). (|2.5p is indeed the starting point of 
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We point out that the above structures, given by the graphs G^, the conductance fields 
Cn and the initial distributions (j,n, could be random. In order to obtain for example a 
quenched hydrodynamic limit, it is enough to prove ()2.5p for almost all realizations of the 
random environment, given by the graph G„ and the conductance field Cn- 



3. Proof of Theorem 12.11 

We write pnit,x,y) for the probability that the random walk Xn{t\x) is in y at time t. 
By the symmetry of the jump rates it holds 

Pn(t,x,y) = pn{t,y,x) , yx,yeVn. (3.1) 

Recall the graphical construction of the simple exclusion process given in the previous 
section and set N^^y for A^"^ ^| . Since 

dr]At)= E {Vy-r].)it-)dN:^yit) 

we can write 

dr]{t) = L„r/(t)(it + dWit), (3.2) 

where 

dM;(t)= iVy-V.){t-)dAly{t), Aly{t) = Nly{t)-c{x,y)t. 

Note that M"(-) has trajectories of bounded variation on finite intervals a.s. 
Formally, (|3.2|) implies that 

rj{t) = T{t)ri{0) + [ T{t- s)dM'^{s) (3.3) 







where T{t) = e*'^", i.e. 



r,^{t)=Y,Pn{t,x,y)riy{Q)+Y, [ p{t - s,x,y)dM^{s). (3.4) 

Due to the graphical construction of the dynamics, if J^j-gy ^^^(O) < ^i then for all but 
a finite family of indexes y dMy{s) = for all ^ s ^ t and in particular the last series 
in ()3.4p reduces to a finite sum and is meaningful. In this case, one can check that ()3.4p 
holds a.s. by direct computation using that 

^Pn{t, X, y) = (L„p„(i, •, y))^ . 

The following result shows that the site exclusion constraint is negligible from a hydrody- 
namic viewpoint: 

Proposition 3.1. Given (5 > 0, t > 0, (/9 E Cc{X) and given a sequence of probability 
measures /u„ on {0, 1}^", it holds 



ntoo 



- Y 'P{x)r]^{t) - — X] ^^^^ X] Pn{t,x,y)r]y{0) > 6) = 0. (3.5) 



xev„ xeVn y&Vn 
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Proof. Let the support of be included in A and fix e > 0. Given x £Vn and t > define 
C^{t) as the connected component of Gt containing x. Then for each positive integer n 
we can choose a bounded set Bn C X such that P{An) < e where An is the subset of 
configurations uj satisfying the following conditions: 



Uxey„nAC(t) CS„, (3.6) 
^ E E Pnit,x,y)^6/2. (3.7) 

x(iV,^n/\y(iVn\B^ 

Given r/(0) and n, we define r/(")(0) G {0,1}^ as rjx '{Q) = rjx{Q)Ix^y^r\B„ and write r]^'^'{s) 
for the configuration at time s obtained by the graphical construction when starting from 
r/(")(0) at time 0. 

Due to the graphical construction of the dynamics and condition (|3.6p , ii uj G An then 
Vxit) = iTx'\'t) for all x G n Bn, thus implying 



;^^9^(x)r,.(t) = ;^^V'(x)4")(t). 



Moreover, due to ()3.7p . if a; € .4^ then 

- E E Pn(t,x,y)|7?,(0) - r?(")(0)I ^ 5/2. 



Therefore the l.h.s. of ()3.5p with fixed n can be bounded by 

{\Zn\> 5/2) + P{Al) ^ 4E^„ {Zl) /5^+e, 

where 

Zn = -Y. ^{x)vt\t) ^(^) E Pn{t,x,yU^\Q). 



Since Y^^&v^ 'ni"'\o) < oo, setting 
p.4p implies that 

= — E '^(^) E / Pn(t-5,X,y)dM^")(s). 

In order to conclude the proof it is enough to apply Lemma 13.21 below to the above 
estimates. □ 

Lemma 3.2. For each n ^ 1 let f„ be a probability measure on {0, 1}^" such that 

i^nCY Vx < oo) = I . 

xeVn 

Then 

1 ^ ^ f' \2 



lim 

ntoo 



[ Pn{t-s,x,y)dM;^{s] 



a,, 

XGVn y&Vn 



0. 



Recall that the above series over y reduces to a finite sum whenever "^^^y^ Vx{0) < oo, 
and therefore it is well defined a.s. 
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Proof. We define fn as 

1 /■* 



x£Vn y&Vn 

2a~ ^ 



1 f 

— 5Z XI X] {Pn{t-s,x,y)-pn{t-s,x,z)){ri^-r]y){s-)dA'^,^{s). 

(3.8) 



We remark that due to the graphical representation of the exclusion process, /„ can be 
thought of as a function on the probability space ({0, 1}^" x il, Bn x J", z/„ (g) P) , where Bn 
denotes the Borel u-algebra of the product space {0, 1}^". Moreover, note that ^ c{(p) 
due to (I33D. 

In the following arguments n can be thought of as fixed. Due to our assumption on f„, 
given £ with < e < 1 there exists a bounded set C„ C X such that ^ e where 

A = {ri : r]^ = if X e Vn\Cn} . 
Moreover, one can find a bounded set D„ C X such that F{B'^) ^ e, where 

Then (z/„ (g) P)(j4 x 5) ^ (1 — e)^. Due to the graphical representation, one gets IaxbIu = 
Iaxb/, where 

/ = 

^ X] '^(^^ Y 12 I {pn{t-s,x,y)-pn{t-s,x,z)){rj^-r]y){s-)dA'^,^{s). 
In particular, 

IE..(/n) < c{^? iyn X P) ((^ X Bf) + (lAxe/') ^ c((^)2(2e - e^) + . 



Since the quadratic variation of the martingale is c(6)t and since all series are finite 
sums, we can compute 



/ ds(r/2 - %)^(s) [ ^ (/^(x) - s, X, y) - - s, x, z)) 



^ E E ^(y^^) ds^J2 fi^)[Pn{t - s,x,y) - pnit - s,x,z) 



^ E E c{y,z)['ds{P^^{y)-P^^{z)f^^fds<P:^,-I.P, 



4^2 



1 /■* d 1 

— / ds— < P^ip^PsUip >^„= — (< V3,<y9 >rr»„ - < Pt^^P't^ >m„) ^ 

2an Jo ds 2an 
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where < •, ■ >m„ denotes the scalar product in L^(y„,m„) (recall ()2.ip ). Note that in the 
second identity we have used the symmetry ()3.ip which implies: 



^ (p{x)(pn{t - S,X,y) -Pnit - S,X,z) \ = Pt-^s'^iv) " Pt-~s'^{z) . 



X&Vn 



The above bound, the vague convergence (j2.ip and the assumption a„ ^ oo imply that 
lim„^oo '^vn [fn] ^ c(99)^(2e — e^). Since e is arbitrary, we get the thesis. □ 

3.1. Conclusion. We have now all the tools to prove Theorem l2.ll Indeed, by Proposition 
3.11 in order to derive the thesis we only need to show that, given 6 > 0, 



lim P. 

nfoo 



- J2 XI Pn{i,x,y)r]y{0) - f ip{u){Ptpo){u)m{du) 



x£Vn y£Vn 

Due to the symmetry ()3.ip we can write 



>6)=0. (3.9) 



1 



Vix) Pn{t,x,y)riy{0) = — Yl Vx{0)Pr^ix) 



xeVn y€Vn XGVn 

By the homogenization assumption (|2.5p . we only need to prove that 



lim fin 

nfoo 



- r]^{0)PMx) - [ ip{u){Ptpo)iu)midu) > 6/2 



(3.10) 



Since W has symmetric kernel, we conclude that 



(p{u){Ptpo){u)m{du) 
Hence, we only need to prove that 



Ptip{u)po{u)m{du) 



hm p,n 



This follows from assumption ()2.6p . 



— V VxiO)PMx) - [ poiu)PMu)m{du) > 5/2) =0. (3.11) 

X&Vn 
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